Universal Dynamical Steps in the Exact Time-Dependent Exchange-Correlation Potential 
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We show that the exact exchange-correlation potential of time-dependent density-functional the- 
ory displays dynamical step structures that have a spatially non-local and time non-local dependence 
on the density. Using one-dimensional two-electron model systems, we illustrate these steps for a 
range of non-equilibrium dynamical situations relevant for modeling of photo-chemical/physical 
processes: field-free evolution of a non-stationary state, resonant local excitation, resonant complete 
charge-transfer, and evolution under an arbitrary field. Lack of these steps in usual approximations 
yield inaccurate dynamics, for example predicting faster dynamics and incomplete charge transfer. 



The vast majority of applications of time-dependent 
density functional theory (TDDFT) today deal with the 
calculation of the linear electronic spectra and response 
of molecules and solids, and provide an unprecedented 
balance between accuracy and efficiency [1, 2]. The the- 
orems of TDDFT also apply to any real-time electron 
dynamics, not necessarily starting in a ground-state, 
and possibly subject to strong or weak time-dependent 
fields. Time-resolved dynamics are particularly impor- 
tant and topical for TDDFT for two reasons. First, there 
is really no alternative practical method for accurately 
describing correlated electron dynamics, and second, 
many fascinating new phenomena and technological 
applications lie in this realm. These include: attosecond 
control of electron dynamics [3], photo-induced coupled 
electron-ion dynamics (for example in describing light- 
harvesting and artificial photosyntheses), and photo- 
chemical/physical processes [4, 5] in general. TDDFT 
in theory yields all observables exactly, solely in terms 
of the time-dependent density however in practice, ap- 
proximations must be made both for the observable 
as a functional of the density, and for the exchange- 
correlation (xc) functional. Thus the question arises as 
to whether the approximate functionals that have been 
successful for excitations predict equally well the dy- 
namics in the more general time-dependent context. In 
particular, the exact xc contribution to the Kohn-Sham 
(KS) potential at time t functionally depends on the his- 
tory of the density n(r, t' < t), the initial interacting 
many-body state "Iv and the choice of the initial KS state 
$o : v xc [n; ^o, $o]( r ; t). However, almost all calculations 
today use an adiabatic approximation, = v^c [ n W]/ 
that inputs the instantaneous density into a ground- 
state xc functional [6, 7], completely neglecting both 
the history- and initial-state-dependence. Further, the 
ground-state functional v^c must be approximated; hy- 
brid functionals, that mix in a fraction of exact-exchange 
to an xc functional that otherwise depends locally or 
semi-locally in space on the density, are most popu- 
lar for the spectra of molecules, while the spatially lo- 
cal LDA and semi-local GGA's are most popular for 



solids (see Ref. [2] and references therein). 

Although understanding when such approximations 
are expected to work well or fail has advanced signif- 
icantly in the linear response regime [2], considerably 
less is known about the performance of approximate 
TDDFT for general non-linear time-dependent dynam- 
ics [8-10]. Part of the reason for this is due to the lack 
of exact, or highly accurate, results to compare with. 
Moreover, even in the case where an accurate calcula- 
tion is available, it is very complicated to extract the ex- 
act xc potential (although see Refs. [11, 12] for significant 
progress). Thus, it is critical for the reliability of TDDFT 
for describing fundamental dynamical processes in the 
applications mentioned earlier, to first test available xc 
approximations on systems for which the exact xc po- 
tential can be extracted. One such case is that of two- 
electrons in a spin-singlet, chosen to start in a KS single- 
Slater determinant. We show that, in this case, the usual 
adiabatic and semi-local approximations typically fail to 
capture a critical and fundamental structure in the ex- 
act correlation potential: a time-dependent step, that has 
a spatially ultranonlocal and non-adiabatic dependence 
on the density. This feature is missing in all available 
TDDFT approximations today. Even the exact adiabatic 
functional misses this dynamical step structure. This 
leads to erroneous dynamics, e.g. faster time scales are 
observed in the adiabatic approximations for examples 
where the step opposes the density evolution. 

For two-electrons in a spin-singlet we choose, as is 
usually done, the initial KS state as a doubly-occupied 
spatial orbital, <f>(r,t). Then the exact KS potential for 
a given density evolution can be found easily [13]. In 
one-dimension (ID), we have 

(d x n(x,t)) 2 d£n(x,t) u 2 (x,t) f x du(x',t) 
Vsi - X ' t >- 8n 2 (x,t) + 4n(x,t) 2 J dt 

(1) 

where u(x, t) = j(x,t)/n(x,t) is the local "velocity", 
n(x, t) is the one-body density, and j(x, t) is the current- 
density. We numerically solve the exact time-dependent 
Schrodinger equation for the two-electron interacting 



wavefunction, obtain n(x, t) and j(x, t), and insert them 
into Eq. 1. The exchange-potential in this case is simply 
minus half the Hartree potential, v x (x, t) = —v H (x, t)/2, 
with v H (x, t) = J w(x' , x)n(x' , t)dx' , in terms of the two- 
particle interaction w(x' , x). Therefore, we can directly 
extract the correlation potential using 

v c (x,i) = v s (x,t) - v eKt (x,t) - v„(x,t)/2 , (2) 

where v CKt {x 7 t) is the external potential applied to the 
system. The two electrons in all our ID examples in- 
teract via the soft-Coulomb interaction [14], w(x',x) = 
1 / \J (x' — x) 2 + 1. We use atomic units throughout. 

We start the analysis with some purely (or largely) 
two-state systems, in which the exact interacting time- 
dependent wavefunction, can be expanded in a 
basis consisting of the ground-state, \^> g (t)), and the first 
excited singlet state, \^ e (t)}: 

|*(f)> =O s (f)|tt fl >+ 0e (t)|* e ) (3) 

where a g (t) and a e (t) are coefficients given by: 
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where E g , E e are the energy eigenvalues of the two 
states, d ab = J **(xi, x 2 )(xi + x 2 ) 1 $'b(xi,X2)dxidx2 is 
the transition dipole moment and £{t) = Acos(ojt) is 
an applied electric field of strength A and frequency tu. 
In the weak amplitude limit, with u 3> |d e9 A| and u 
close to the resonant frequency, this reduces to the text- 
book Rabi problem. When on-resonance, the system os- 
cillates from one state to the other over a Rabi cycle of 
period Tr = 2tt / (\d eg \A) . By solving Eq. 4 we can eas- 
ily construct the current and density at any time, their 
time-derivatives, and hence all pieces entering Eq. 1. 

In our first example, we consider a "ID He atom", 
where 7» cx t = — 2j\Jx 2 + 1, subject to a weak electric field 
of strength A = 0.00667au and frequency u = 0.533au, 
resonant with the first singlet excitation [15, 16]. Fig- 
ure 1 shows snapshots of the correlation potential over 
one Rabi cycle, while Fig. 2 shows snapshots over one 
optical period centered around Tr/4. (Note that the sys- 
tem is not exactly periodic over the Rabi period as the 
two time-scales dictated by the optical frequency and 
the Rabi frequency are not commensurate). 

The most salient feature of the correlation potential is 
the presence of time-dependent steps, that oscillate on 
the time-scale of the optical field. These steps arise from 
the fourth term of Eq. 1: whenever there is a net "acceler- 
ation", dtu(x, t), through the system, the spatial-integral 
is finite, resulting in a potential rising from one end of 
the system to the other. The correlation potential thus 
has a spatially ultranonlocal dependence on the density, 
as it changes far from the system. 

Further, the time-dependence of the steps is non- 
adiabatic, meaning that the instantaneous density is 
not enough to determine the correlation potential func- 
tional. This is clear from Figure 2 where the small 
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FIG. 1: (color online). Snapshots of the exact correlation po- 
tential (solid black), density (blue dotted), and exact-adiabatic 
(red dashed) over one Rabi cycle; at T R /2 the density of the 
first excited state is essentially exactly reached. In all graphs 
in this paper, the correlation potentials are plotted up to an 
arbitrary irrelevant time-dependent constant. 




FIG. 2: (color online) Snapshots of the correlation potential 
(left), and corresponding density (right), at times indicated in 
the right panel. 



changes of the density between time steps cannot cap- 
ture the observed large changes of the step. One is 
tempted to point to the time-derivatives in the fourth 
term in Eq. 1 as further evidence for the non-adiabatic 
dependence, however caution would be needed for such 
an argument as time non-locality in v s is not the same 
as time non-locality in v c [2]: the fourth term, for ex- 
ample, may be written as v cxt plus other terms, and 
although v cxt has typically strongly non-adiabatic de- 
pendence, this is irrelevant because it is never approx- 
imated as a functional in practice [2, 17], rather it is 
taken from the problem at hand. Only the xc potential 
must be approximated, and the functional-dependence 
of this cannot be deduced directly from Eq. (1). In- 
stead, to unambiguously show the non-adiabatic de- 
pendence of the step, we plot the "adiabatically-exact" 
correlation potential in Fig. 1. This is defined by the 
exact correlation potential for which both the interact- 
ing and KS wavefunctions are ground-state wavefunc- 
tions with density equal to the instantaneous one i.e. 
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v adia-ex^ = V »di»[ n ]- V ^[ n ]- Vu [ n ]- Vx [ n ] [18], where 
Vg * a [ n ] is the external potential for two interacting elec- 
trons whose ground-state has density n, and v^ dia [n] 
is the exact ground-state KS potential for this density 
(given by the first two terms in Eq.(l). We find v^l a [n] 
using similar techniques to Ref. 18 (see also Ref. 19 ). 
Fig. 1 shows that v ^ dla - ex [j-j] indeed does not capture the 
dynamical step structure. 

Before turning to our next example, we verify that 
the two-state approximation for the dynamics is accu- 
rate enough for our purposes. Actually one aspect of 
the potentials we find is indeed an artifact of the two- 
state approximation: the correlation potential asymp- 
totically has a slope so to exactly cancel the externally 
applied electric field. This is because the two-state ap- 
proximation cannot correctly describe polarization aris- 
ing from occupying many excited states in time. The 
KS potential obtained from inverting the two-state ap- 
proximation must therefore be flat asymptotically, as it 
cannot describe states that are polarized asymptotically. 
The field is so weak in our case that this effect is hardly 
noticeable on the scale of Figs. 1 and 2, but to check 
that our conclusions regarding the dynamics step struc- 
ture are unaffected by the two-level approximation, we 
computed the KS potential using the density, current, 
and their time-derivatives from the numerically exact 
wavefunction, found using octopus [20, 21]. Apart 
from some extra structure in the tail region (small peaks 
and steps as we move away from the atom), and the 
linear field-counteracting term, the correlation potential 
agrees with that from the two-state model. 

Dynamical step features have arisen in TDDFT in ear- 
lier studies; Refs. [18, 22] showed they appear in ion- 
ization processes, and linked them to a time-dependent 
derivative discontinuity, related to fractional charges. 
In time-resolved transport, step structures have been 
shown to be essential for describing Coulomb-blockade 
phenomena [23], again related to the discontinuity. In 
the response regime, field-counteracting steps develop 
across long-range molecules [24]. In open-systems- 
TDDFT, Ref. [25] shows steps arise when using a closed 
KS system to model an open interacting one. In the lin- 
ear response regime, the xc kernel for charge-transfer 
excitations displays frequency-dependent steps [26]. 
However we argue that the dynamical step structures 
we are seeing are generic, and moreover, unlike most of 
the above cases [18, 22-24], cannot be captured by an 
adiabatic approximation. They appear with no need for 
ionization nor subsystems of fractional charge, nor any 
applied field (see next example), unlike in Refs. [18, 22- 
24]. In this sense our results are more akin to Ref. [11], 
which studies the physically very different situation 
when an electron freely propagates through a wire. The 
range of the examples we present suggests that such 
non-adiabatic and non-local steps generically arise when 
dealing with real electron dynamics. 

Our second example accentuates the fact that dynam- 
ical step structures need neither ionization nor an exter- 
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FIG. 3: (color online) The exact correlation potential (black 
solid) shown at times 0, 2, 4, and 6 au, for the two-state ex- 
ample (Eq. 5). Also shown are the local acceleration (red dot- 
ted), f x dtu — t'cxt (purple dashed), and the adiabatically-exact 
correlation potential (blue dash-dot). 



nal field to appear. We begin in an equal linear superpo- 
sition of the ground and first-excited state of the ID He 
and let it evolve freely, so that 

|*(t)) = (e-^'l^) + e- iE ^ e )) /V2 . (5) 

It will oscillate back and forth between the two states 
with frequency luq = E e — E g . The two-state approxi- 
mation for this is exact at all times. Again, we see large 
steps in the correlation potential, as shown in Fig. 3. 
To support the discussion and provide a microscopic in- 
sight behind this phenomenon, we also plot in Fig. 3 the 
acceleration, a(x,t) = d t u(x,t), and its spatial integral 
with the external potential subtracted out. The position 
and magnitude of the step at each time is heavily de- 
pendent on this term. Peaks in the acceleration, when 
integrated, become local steps in the potential and the 
asymptotic value of the step in v c (x, t) is given by the 
total step in the spatial integral of a(x, t). Although lo- 
cal step-like features may be cancelled out by the other 
terms in Eq. (1), the net magnitude of the step is deter- 
mined from the asymptotic values of this integral. 

Note that we have the freedom to choose the initial 
state of the KS system as long as it has the same density 
and first derivative in time of the exact density [27], and 
the shape of the exact correlation potential depends on 
this choice [19]. We used a doubly occupied orbital in 
the previous example, so that we can calculate the exact 
v c using the method discussed. A different choice, with 
a configuration more similar to that of the interacting 
initial state could well yield a more gentle correlation 
potential [19], with less dramatic step structure. 

The generality of the step feature in dynamics is fur- 
ther supported by considering different resonant excita- 



3 



+2Topl/3 



FIG. 4: (color online) The correlation potential (left) and den- 
sity (right) shown at snapshots of time Tr/8± fractions of the 
optical period, T op t = 27r/0.112au, that are indicated in the 
key, for the two-well potential model Eq. (6). 
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tions. Consider a double-well as a model of a molecule: 
2 1 

VeKt(X,t) = 7£ + 3.5)2 + Lu " C0S h 2 (x - 3.5) ~ £{t)X 

(6) 

with £(t) = 0.006 cos(O.im). Here the ground-state has 
two electrons in one well, and a charge-transfer excited 
state f e , with one electron in each well, at a frequency 
of 0.112au. We use the ground-state and ^> e in the two- 
state model of Eq. (3), and solve for the occupations us- 
ing Eq. (4); we again checked the two-state result against 
the exact numerical solution using octopus. The system 
behaves like the Rabi problem with non-zero dipole mo- 
ment for the ground-state [28, 29]. In Fig. 4 we plot the 
correlation potential for several times within an optical 
cycle around Tr/8. Again dynamical steps oscillating 
on the optical frequency time scale emerge. The situa- 
tion is more complicated as a step related to the der- 
ealization of the density during the charge-transfer pro- 
cess slowly develops (on the time-scale of Tr/2) [30]. 
The dynamical step can then increase, decrease, or even 
reverse this charge-transfer step. Approximations un- 
able to develop steps lead to incomplete charge-transfer. 
This, along with other details of time resolved charge- 
transfer, is investigated in more detail in Ref. [30]. Eor 
our purposes, it is sufficient to note that dynamical steps 
are again present to capture the exact dynamics. 

Finally, we explicitly demonstrate that the non-local 
non-adiabatic step feature is a generic aspect of the cor- 
relation potential in the following way. We subject the 
ID He atom to an electric field that is chosen somewhat 
arbitrarily: it is relatively strong and linearly switched 
on over 2 optical cycles, with an off -resonant frequency. 
In Fig. 5, we show the exact correlation potential at 
four times, along with the density, and the adiabatically- 
exact correlation potential. The time-dependent step in 
the exact v c is once again evident, and again fails to be 
captured by the adiabatically-exact approximation. 

In summary, dynamical steps in the correlation poten- 
tial are a generic feature of electron dynamics. The step 
features discussed above arise from part of the fourth 
term of Eq. (1), which suggests that any time there is 
a net localized acceleration, there will be a step, and 



FIG. 5: (color online). The exact correlation potentials (solid, 
black) at times 10,15,20,25au during propagation under £{t) = 
^775 ^jf cos(0.3 t). Also shown are the adiabatically-exact cor- 
relation potentials (dashed red), and the density (dotted blue, 
scale on the right). 



that it will have a very non-local spatial dependence 
on the density, and is non-adiabatic. This represents a 
type of time-dependent screening, where the electron- 
electron interaction hinders electron movement to cer- 
tain regions. Although two-electron systems were stud- 
ied here, we expect steps are a more general feature of 
electron dynamics, as supported by the recent Ref. [11]. 

The lack of the step in approximations leads to in- 
correct dynamics. For example, faster time scales in 
adiabatic approximations were found for the field-free 
dynamics of a linear superposition state, where the di- 
rection of the step tended to oppose the density's mo- 
tion. The exact dipole and adiabatic exact-exchange 
(AEXX) dipole for this case are shown in Fig. 6. We com- 
puted the dynamics of the the local excitation in Figs 1 
and 2 using AEXX, adiabatic LDA, and adiabatic self- 
interaction-corrected LDA. In all cases, we found the 
timescale for the dipole oscillations was faster than in 
the true case, in spite of providing good linear response 
spectra [16]. How general this finding is will be investi- 
gated in future work. 

We note that the xc electric field, defined as the gra- 
dient of the xc potential, has a more local character than 
the potential. This suggests that considering functional 
approximations to this field, or, more generally, to an xc 
vector potential [11, 31], may point to an easier path to 
develop approximations containing step features. 

As applications of TDDFT continue to expand, it is 
crucial to further study what the impact of the miss- 
ing steps in the approximations are on the predictions 
of these calculations. When starting in the ground-state, 
the exact adiabatic potential may follow well the exact 
dynamics at short times, but as soon as there is an ap- 
preciable change in the occupation of an excited state, 
the exact soution develops the dynamical step, entirely 
missing in the adiabatic one. This result is general and 
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FIG. 6: (color online). The dipole moments in the field-free 
propagation of the linear-superposition state (see text), where 
the TDKS calculation starts in the exact singlet doubly occu- 
pied orbital initial state. 

applies to all available functionals. It raises an impor- 
tant issue when applying TDDFT to fundamental photo- 
induced processes (e.g. photovoltaics, artificial photo- 



synthesis, photoactivated chemistry, photophysics, etc): 
all these involve a significant change of state population. 
Clearly population of many-body states due to the exter- 
nal field is not a linear process and requires functionals 
able to cope with the generic features of the dynamical 
step that we have unveiled in the present work. 
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